MONTE CARLO STUDY OF SUPERNOVA NEUTRINO 

SPECTRA FORMATION 



o 
o 

(N 
X> 

IX, 



Mathias Th. Keil and Georg G. Raffelt 

Max-Planck-Institut fur Physik ( Werner- Heisenberg-Institut) 
Fdhringer Ring 6, 80805 Miinchen, Germany 

Hans-Thomas Janka 

Max-Planck-Institut fur Astrophysik 
Karl-Scharzschild-Str. 1, 85741 Garching, Germany 



(N 

(N 
> 

in 
m 
o 

oo 
o 

O 

a: 
6 

c3 



ABSTRACT 

The neutrino flux and spectra formation in a supernova core is studied by using a Monte 
Carlo code. The dominant opacity contribution for is elastic scattering on nucleons i^N — > 
Nz^, where always stands for either or v T . In addition we switch on or off a variety of 
processes which allow for the exchange of energy or the creation and destruction of neutrino pairs, 



notably nucleon bremsstrahlung NN 



NNi/ M z7„, 



the pair annihilation processes e + e 



and 



v c v c —> VuVn, recoil and weak magnetism in elastic nucleon scattering, elastic scattering on 



electrons 



and elastic scattering on electron neutrinos and anti- neutrinos v^v c 



v c v^ and v^v^ — > v^v^. The least important processes arc neutrino- neutrino scattering and 
e + e _ annihilation. The formation of the spectra and fluxes of is dominated by the nucleonic 
processes, i.e. bremsstrahlung and elastic scattering with recoil, but also v e i> annihilation and 
scattering contribute significantly. When all processes are included, the spectral shape of 
the emitted neutrino flux is always "pinched," i.e. the width of the spectrum is smaller than that 
of a thermal spectrum with the same average energy. In all of our cases we find that the average 
energy exceeds the average v c energy by only a small amount, 10% being a typical number. 
Weak magnetism effects cause the opacity of to differ slightly from that of , translating into 
differences of the luminosities and average energies of a few percent. Depending on the density, 
temperature, and composition profile, the flavor-dependent luminosities L Ue , L Pa , and L v can 
mutually differ from each other by up to a factor of two in either direction. 

Subject headings: diffusion — neutrinos — supernovae: general 



1. INTRODUCTION 

In numerical core-collapse supernova (SN) sim- 
ulations, the transport of \i- and r-neutrinos has 
received scant attention because their exact fluxes 
and spectra are probably not crucial for the explo- 
sion mechanism. However, the recent experimen- 
tal evidence for neutrino oscillations implies that 
the flavor-dependent fluxes and spectra emitted 
by a SN will be partly swapped so that at any dis- 
tance from the source the actual fluxes and spec- 
tra can be very different from those originally pro- 



duced. In principle, this effect can be important 
for the SN shock revival (Fuller et al. 1992) and 
r-process nucleosynthesis (Qian et al. 1993, Pastor 
& Raffelt 2002), although the experimentally fa- 
vored small neutrino mass differences suggest that 
this is not the case. On the other hand, in view 
of the large- mixing- angle solution of the solar neu- 
trino problem flavor oscillations are quite relevant 
for the interpretation of the SN I987A neutrino 
signal (Jegerlehner, Neubig, & Raffelt 1996, Lu- 
nardini & Smirnov 2001a, Kachelriess et al. 2002, 
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Smirnov, Spergel, & Bahcall 1994). More im- 
portantly, the high-statistics neutrino signal from 
a future galactic SN may allow one to differen- 
tiate between some of the neutrino mixing sce- 
narios which explain the presently available data 
(Chiu & Kuo 2000, Dighe & Smirnov 2000, Dutta 
et al. 2000, Fuller, Haxton, & McLaughlin 1999, 
Lunardini & Smirnov 2001b, 2003, Minakata & 
Nunokawa 2001, Takahashi & Sato 2002). Even 
though the solution of the solar neutrino problem 
has been established, the magnitude of the small 
mixing angle 613 and the question if the neutrino 
mass hierarchy is normal or inverted will remain 
open and can be settled only by future precision 
measurements at dedicated long-baseline oscilla- 
tion experiments (Barger et al. 2001, Cervera et al. 
2000, Freund, Hubcr, & Lindner 2001) and/or the 
observation of a future galactic SN. 

The usefulness of SN neutrinos for diagnosing 
flavor oscillations depends on the flavor depen- 
dence of the fluxes and spectra at the source. Very 
crudely, a SN core is a black-body source of neutri- 
nos of all flavors which are emitted from the sur- 
face of the proto-neutron star that was born after 
collapse. It is the flavor-dependent details of the 
neutrino transport in the neutron-star atmosphere 
which cause the spectral and flux differences that 
can lead to interesting oscillation effects. 

The v c and v c opacity is dominated by the 
charged-current processes c c n — > pe~ and D c j> — > 
ne + , reactions that allow for the exchange of en- 
ergy and lepton number between the medium and 
the neutrinos. Therefore, it is straightforward to 
define an energy-dependent neutrinosphere where 
this reaction freezes out for neutrinos of a particu- 
lar energy. This sphere yields a thermal contribu- 
tion to the neutrino flux at the considered energy. 
The atmosphere of the proto-neutron star is neu- 
tron rich, providing for a larger v c opacity than 
for D c so that for a given energy the v e flux orig- 
inates at deeper and thus hotter layers than the 
v c flux. In other words, a larger fraction of the v c 
flux emerges with high energies. This simple ob- 
servation explains the usual hierarchy (e Po ) > (e„ e ) 
of the mean energies. The spectra are found to 
be "pinched" , meaning that the high-energy tail 
is suppressed relative to that of a thermal spec- 
trum with the same mean energy (Janka & Hillc- 
brandt 1989a, b). This numerical result can be un- 
derstood analytically by constructing the neutrino 



spectrum from the fluxes emitted by the energy- 
dependent neutrinospheres which are at different 
temperatures (Myra, Lattimer, & Yahil 1988, Gio- 
vanoni, Ellison, & Bruenn 1989). 

The formation of the v^, 9^, v T , and 9 T spec- 
tra is far more complicated. The opacity is dom- 
inated by the neutral-current scattering on nucle- 
ons, i^N — > Nz/^, a process that prevents neutrino 
free streaming, but is unable to change the neu- 
trino number and is usually considered to be in- 
efficient at exchanging energy. (Here and in the 
following stands for either or v T .) Neutrino 
pairs can be created by nucleon bremsstrahlung, 
NN — > NNi^P^, and pair annihilation, e~e + — > 
Vpf/p or v c v c — > Vpi/p, while v^v^ pairs are ab- 
sorbed by the inverse reactions. In addition, en- 
ergy is exchanged by elastic scattering on leptons, 
notably i/^e - — > e~v /1 , by the recoil in nucleon 
scattering, i/ M N — > N^, and by inelastic scatter- 
ing on nucleons ^ M NN — ► NN^ M , a channel that 
is the "crossed process" of bremsstrahlung. For 
a given neutrino energy these processes freeze out 
at different radii so that one can define a "number 
sphere" for the pair processes, an "energy sphere" 
for the energy-exchange processes, and a "trans- 
port sphere" for clastic nucleon scattering with 

-Rnumbcr < -Ronorgy < -^transport (Suzuki 1990). 

The region between the number sphere and the 
transport sphere plays the role of a scattering at- 
mosphere because neutrinos can not be created 
or destroyed. They propagate by diffusion and 
can still exchange energy with the background 
medium. 

Usually the transport sphere is deeper than 
the v e sphere so that numerical simulations find 
( £ iv) > ( ep c) > ( e ^e)- This hierarchy is the main 
motivation for the proposed use of SN neutrinos 
cXS (X diagnostic for neutrino oscillations. However, 
the quantitative statements found in the literature 
range from (e„ ) being 20% to nearly a factor of 2 
larger than (ep c ); for a review see Janka (1993) and 
Sec. 4.3. Of course, the mean energies and their 
ratios change significantly between the SN bounce, 
accretion phase, and the later neutron-star cool- 
ing phase. Therefore, one must distinguish care- 
fully between instantaneous fluxes and spectra and 
the time-integrated values. While for the analysis 
of the sparse SN 1987A data only time-integrated 
values make sense, a future galactic SN may well 
produce enough events to study the instantaneous 
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fluxes and spectra (Barger, Marfatia, & Wood 
2001, Minakata et al. 2001). 

The overall energy emitted by a SN is often 
said to be equipartitioned among all six neutrino 
degrees of freedom. In some numerical simulations 
the neutrino luminosities are indeed astonishingly 
equal for all flavors (Totani et al. 1998), while 
other simulations easily find a factor of two dif- 
ference between, say, the P M and v c luminosities, 
at least during the accretion phase (Mezzacappa 
et al. 2001). Therefore, it is by no means obvious 
how precisely equipartition can be assumed for the 
purpose of diagnosing neutrino oscillations. 

Another important feature is the neutrino spec- 
tral shape, notably the amount of pinching. If 
one could assume with confidence that the instan- 
taneous spectra of all flavors arc pinched at the 
source, and if the measured SN neutrino spec- 
tra were instead found to be anti-pinched, this 
effect would be a powerful diagnostic for the par- 
tial spectral swapping caused by flavor oscillations 
(Dighe & Smirnov 2000). 

Unfortunately, the existing literature does not 
allow one to develop a clear view on these "fine 
points" of the neutrino fluxes and spectra, largely 
because not enough attention has been paid to 
the and v T emission from a SN core. The 
published full numerical SN collapse simulations 
have not yet included the bremsstrahlung process 
or nuclcon recoils (but see first results of state- 
of-the-art models in Rampp et al. 2002), even 
though it is no longer controversial that these ef- 
fects are important (Janka et al. 1996, Burrows 
et al. 2000, Hannestad & Raffelt 1998, Raffelt 
2001, Suzuki 1991, 1993, Thompson, Burrows, & 
Horvath 2000). Moreover, some of the interest- 
ing information such as the spectral pinching was 
usually not documented. 

Another problem with self-consistent hydrody- 
namic simulations is that the models with the 
most elaborate neutrino transport usually do not 
explode so that even the most recent state-of- 
the-art simulations do not reach beyond the ac- 
cretion phase at a few hundred milliseconds af- 
ter bounce (Rampp & Janka 2000, Mezzacappa 
et al. 2001, Liebendorfer et al. 2001), thus not pro- 
viding any information on the neutron-star cool- 
ing phase. Successful multi-dimensional models of 
the explosion (e.g., Fryer & Warren 2002, Fryer 
1999 and references therein) were also not con- 



tinued to the neutron-star cooling phase. These 
simulations, moreover, treat the neutrino trans- 
port only in a very approximate way and do not 
provide spectral information. The calculations 
performed by the Livermore group also yield ro- 
bust explosions (Totani et al. 1998). They include 
a mixing-length treatment of the phenomenon 
of neutron-finger convection in the neutron star, 
that increases the early neutrino luminosities and 
thus enhances the energy transfer by neutrinos to 
the postshock medium (Wilson & Mayle 1993). 
Whether neutron-finger convection actually oc- 
curs inside the neutrinosphcre and has effects on a 
macroscopic scale, however, is an unsettled issue. 

We will follow here an alternative approach to 
full hydrodynamic simulations, i.e. we will study 
neutrino transport on the background of an as- 
sumed neutron-star atmosphere. While this ap- 
proach lacks hydrodynamic self-consistency, it has 
the great advantage of allowing one to study sys- 
tematically the influence of various pieces of mi- 
croscopic input physics and of the medium pro- 
file. The goal is to develop a clearer picture of the 
generic properties of the SN neutrino spectra and 
fluxes and what they depend upon. 

To this end we have adapted the Monte Carlo 
code of Janka (1987, 1991) and added new mi- 
crophysics to it. We go beyond the work of 
Janka & Hillebrandt (1989a, b) in that we in- 
clude the bremsstrahlung process, nucleon re- 
coils and weak magnetism, v c D c pair annihila- 
tion into Vfii>^, and scattering of on v e and D e . 
With these extensions we investigate the neutrino 
transport systematically for a variety of medium 
profiles that are representative for different SN 
phases. One of us (Raffelt 2001) has recently 
studied the spectra-formation problem with 
the limitation to nucleonic processes (elastic and 
inelastic scattering, recoils, bremsstrahlung), to 
Maxwcll-Boltzmann statistics for the neutrinos, 
and plane-parallel geometry. Our present study 
complements this more schematic work by includ- 
ing the leptonic processes, Fermi-Dirac statistics, 
and spherical geometry. In addition we apply our 
Monte Carlo code to the transport of v c and v e 
and thus are able to compare the flavor-dependent 
fluxes and spectra. 

In Sec. 2 we first assess the relative importance 
of different processes in terms of their energy- 
dependent "thermalization depth". In this con- 
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text we introduce a number of stellar background 
models. In Sec. 3 we perform a Monte Carlo study 
of transport on the previously introduced back- 
ground models in order to assess the importance 
of different pieces of input physics. In Sec. 4 we 
compare the fluxes and spectra with those of 
v e and v e . We conclude in Sec. 5 with a discussion 
and summary of our findings. 

2. THERMALIZATION DEPTH OF 

ENERGY-EXCHANGE PROCESSES 

2.1. Simple Picture of Spectra Formation 

One of our goals is to assess the relative im- 
portance of different neutrino interaction chan- 
nels with the background medium of the SN core. 
As a first step it is instructive to study the ther- 
malization depth of various energy-exchange pro- 
cesses. Within the transport sphere, the neutri- 
nos are trapped by elastic scatterings on nucleons, 
i^N — > N^ M , which are by far the most frequent 
reactions between neutrinos and particles of the 
stellar medium. (Unless otherwise noted "neutri- 
nos" always refers to any of v^, v^, v T or v T .) 
Assuming for the moment that these collisions are 
iso-energetic (no nucleon recoils), it is straightfor- 
ward to define for a neutrino of given energy e 
the location ( "thcrmalization depth" ) where it last 
exchanged energy with the medium by a reac- 
tion such as v^e~ — ► e~u^. Following Shapiro & 
Tcukolsky (1983) we define the optical depth for 
energy exchange or thermalization by 



Tthcrm(r) 



dr'i 



1 



\ E (r') 



1 



1 



Mr') X E (r') 



(1) 

Here, A_e is the mean free path (mfp) for the rel- 
evant energy-exchange process and At the trans- 
port mfp, i.e. the mfp corresponding to the cross 
section for momentum exchange in the f N — ► Nzv 
reaction. The quantities Ttherm, and At are all 
understood to depend on the neutrino energy e. 
The main philosophy of Eq. (1) is that a neutrino 
trapped by elastic scattering has a chance to ex- 
change energy corresponding to its actual diffu- 
sive path through the scattering atmosphere; for 
a discussion see Suzuki (1990). The thermaliza- 
tion depth -Rthorm is given by 



Ttherm (-^therm) 



2 

3 ' 



(2) 



where i?thcrm depends on the neutrino energy e. 

When this energy dependence is not too steep 
it makes sense to define an average thermalization 
depth, i.e. an "energy sphere" that for pair creat- 
ing processes is equal to the "number sphere." For 
nucleon bremsstrahlung this requirement is well 
fulfilled (Raffelt 2001) so that one may picture 
the energy sphere as a blackbody surface that in- 
jects neutrinos into the scattering atmosphere and 
absorbs those scattered back. The neutrino flux 
and spectrum emerging from the transport sphere 
is then easily understood in terms of the energy- 
dependent transmission probability of the black- 
body spectrum launched at the energy sphere. 
The transport cross section scales as e 2 , implying 
that the transmitted flux spectrum is shifted to 
lower energies relative to the temperature at the 
energy sphere. This simple "filter effect" accounts 
surprisingly well for the emerging flux spectrum 
(Raffelt 2001). For typical conditions the mean 
flux energies are 50-60% of those corresponding 
to the blackbody conditions at the energy sphere. 

Moreover, it is straightforward to understand 
that the effective temperature of the emerging flux 
spectrum is not overly sensitive to the exact loca- 
tion of the energy sphere. If the energy-exchange 
reaction is somewhat more effective, the energy 
sphere is at a larger radius with a lower medium 
temperature. However, the scattering atmosphere 
has a smaller optical depth so that the higher- 
energy neutrinos are less suppressed by the fil- 
ter effect, partly compensating the smaller energy- 
sphere temperature. For typical situations Raffelt 
(2001) found that changing the bremsstrahlung 
rate by a factor of 3 would change the emerging 
neutrino energies only by some 10%. This finding 
suggests that the emitted average neutrino energy 
is not overly sensitive to the details of the energy- 
exchange processes. 

2.2. Neutron-Star Atmospheres 

In order to determine the location of the ther- 
malization depth for different processes we need 
to define our assumed neutron-star atmospheres. 
As a first example we use a model taken from a 
full hydrodynamic simulation. This model is rep- 
resentative for the accretion phase; henceforth we 
will refer to it as the "Accretion-Phase Model I" 
(Fig. 1). It was provided to us by O. E. B. Messer 
and was already used in Raffelt (2001) for a more 
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Table 1 

Characteristics of power-law models. 





Steep 


Shallow 


p 


10 


5 


q 


2.5 


1 


q/p 


0.25 


0.2 


po [10 14 gem- 3 ] 


2.0 


0.2 


T [MeV] 


31.66 


20.0 


r [km] 


10 


10 



schematic study. Based on the Woosley & Weaver 
15 Mq progenitor model labeled sl5s7b, the New- 
tonian collapse simulation was performed with the 
SN code developed by Mezzacappa et al. (2001). 
The snapshot is taken at 324 ms after bounce when 
the shock is at about 120 km, i.e. the star still ac- 
cretes matter. In this simulation the traditional 
microphysics for transport was included, i.e. 
iso-energetic scattering on nucleons, e+e^ annihi- 



lation and v^e scattering. 

As another self-consistent example (Accretion- 
Phase Model II) we obtained a 150 ms postbounce 
model from M. Rampp (personal communication) 
that uses a very similar progenitor (sl5s7b2). The 
simulation includes an approximate general rel- 
ativistic treatment in spherical symmetry as de- 
scribed by Rampp & Janka (2002). The three 
neutrino flavors are transported with all relevant 
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Fig. 1. — Accretion-Phase Model I, a SN model 
324 ms after bounce from a Newtonian calculation 
(O.E.B. Mcsscr, personal communication). 



Fig. 2. — Accretion-Phase Model II, a SN core at 
150 ms postbounce from a general-relativistic sim- 
ulation. (M. Rampp, personal communication). 
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interactions except v c v c pair annihilation to v^v^ 
(see also Sec. 4.1 and Rampp et al. 2002). 

As another set of examples we use two power- 
law profiles of the form 

T = T„^)\ (3, 

with a constant electron fraction per baryon Y e . 
We adjust parameters such that (e) w 20-25 MeV 
for the emerging neutrinos to obtain model at- 
mospheres in the ballpark of results from proto- 
neutron star evolution calculations. We define 
a "steep" power-law model, corresponding to the 
one used by Raffelt (2001), and a "shallow" one; 
the characteristics are given in Table 1. The shal- 
low model could be characteristic of a SN core 
during the accretion phase while the steep model 
is more characteristic for the neutron-star cooling 
phase. The constant electron fraction Y c is an- 
other parameter that allows us to investigate the 
relative importance of the leptonic processes as a 
function of the assumed Y c . 

2.3. Thermalization Depth 

We now calculate the thermalization depth as 
a function of the neutrino energy e for several 
energy-exchanging processes and the neutron-star 
atmospheres described above. We consider the 
neutrino mfp for nucleon bremsstrahlung NN — > 
NN^f-p, pair annihilation c + c~ — ► v^v^ and 
v c v c — ► v^D^, and scattering on charged leptons 
v^e^ — > e ± f /i . The numerical implementation of 
the reaction rates is described in Appendix B. 

In Figs. 3 and 4 we give the thermalization 
depth -Rthcrm as a function of neutrino energy 
e for the two hydrodynamically self-consistent 
accretion-phase models. From top to bottom the 
panels show the results for v e , D c , and v^, respec- 
tively. The step-like curves represent the tempera- 
ture profiles in terms of the mean neutrino energy, 
(e„) = 3.15 T for non-degenerate neutrinos at the 
local medium temperature; the steps correspond 
to the radial zones of our Monte Carlo simulation. 
The other curves represent i?therm for bremsstrah- 
lung (b), e + e~ annihilation (p), v c D c annihilation 
(n), and scattering on e^ 1 (s). In the case of v c 
and v e we do not include bremsstrahlung and v e v e 
annihilation. Particle creation is dominated by 
the charged current reactions on nucleons (urea). 




20 30 40 50 60 70 

R therm [ km ] 

Fig. 3. — i?thorm as a function of neutrino energy 
e for our Accretion-Phase Model I. From top to 
bottom the panels show the results for v e , v e , and 
Vp. Energy exchanging processes: bremsstrahlung 
(solid line), e + c~ annihilation (dashed), v c D c an- 
nihilation (dotted), and scattering on (dash- 
dotted). "Urea" denotes the charged-current reac- 
tion of v e and D e on nucleons. The steps represent 
(e) = 3.15 T. 

For the power-law models we show i?thcrm for ^ 
in Figs. 5 and 6. The different panels correspond 
to the indicated values of the electron fraction Y e . 
Note that Y e represents the net electron density 
per baryon, i.e. the e _ density minus that of e + so 
that Y e = implies that there is an equal thermal 
population of e~ and e + . 

The absorption rate for the bremsstrah- 
lung process varies approximately as e -1 , the 
transport cross section as e 2 so that the inverse 
mfp for thermalization varies only as e 1 / 2 . This 
explains why -Rthcrm for bremsstrahlung is indeed 
quite independent of e. Therefore, bremsstrahlung 
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Fig. 4. — Same as Fig. 3 for the Accretion-Phase 
Model II 



alone allows one to specify a rather well-defined 
energy sphere. The other processes depend much 
more sensitively on e so that a mean energy sphere 
is much less well defined. 

Both electron scattering and the leptonic pair 
processes are so ineffective at low energies that 
true local thermodynamic equilibrium (LTE) can 
not be established even for astonishingly deep lo- 
cations. Brcmsstrahlung easily "plugs" this low- 
energy hole so that one can indeed expect LTE for 
all relevant neutrino energies below a certain ra- 
dius. For higher energies, the leptonic processes 
dominate and shift the energy sphere to larger 
radii than bremsstrahlung alone. The relative im- 
portance of the various processes depends on the 
density and temperature profiles as well as Y e . 

To assess the role of the various processes for 
the overall spectra formation one needs to spec- 
ify some typical neutrino energy. One possibility 
would be (e) for neutrinos in LTE. Another pos- 




Y e = 0.00 



Y e = 0.05 



Y e = 0.50 



12 13 14 15 16 17 18 19 20 

R therm t km ] 

Fig. 5. — -Rtherm for in the steep power-law 
model with the indicated values of Y c . This figure 
corresponds to the bottom panel of Fig. 3. 

sibility is the mean energy of the neutrino flux, 
in particular the mean energy of those neutrinos 
which actually leave the star. For our power-law 
atmospheres this is always around 20-25 MeV. 
Therefore, the process with the largest i?therm in 
this energy band is the one most relevant for de- 
termining the emerging neutrino spectrum. It ap- 
pears that at least for steep profiles pair annihi- 
lation is never crucial once bremsstrahlung is in- 
cluded, i.e. we would guess that including pair 
annihilation will not affect the emerging neutrino 
spectra. The relevance of electron scattering is 
far more difficult to guess. On the one hand it 
surely is more important than recoil in nucleon 
scatterings for some of the relevant energies, on 
the other hand we are not able to define an energy 
sphere for nucleon recoils because this process is 
different from the others in that neutrinos transfer 
only a small fraction of their energy per scatter- 
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Fig. 6. — Same as Fig. 5 for the shallow power-law 
model. 



ing. Therefore, it is not straightforward to assess 
the relevance of electron scattering compared with 
nucleon recoils on the basis of the various thermal- 
ization spheres alone. 

3. MONTE CARLO STUDY OF MUON 
NEUTRINO TRANSPORT 

3.1. Spectral Characteristics 

In order to characterize the neutrino spectra 
and fluxes emerging from a neutron star we need 
to introduce some simple and intuitive parameters. 
One is the mean energy 



Jo°° dee I-i d ^f( e ^) 



(e) = 



Jo°° de S-i d M/(e,M) 



(4) 



where f{e,fi) is the neutrino distribution function 
with e the energy and \i the cosine of the angle 
between the neutrino momentum and the radial 



direction. If the neutrinos are in LTE without a 
chemical potential one has 



/(e,M) 



1 + cxp(e/T) 



and therefore 



M T 



3.1514T. 



(5) 



(6) 



One can define an effective neutrino temperature 
for non-equilibrium distributions by inverting this 
relationship. 

It is often useful to extract spectral characteris- 
tics for those neutrinos which are actually flowing 
by removing the isotropic part of the distribution. 
Specifically, we define the average flux energy by 



io°° de /-i'^M/^m) 



(7) 



Far away from the star all neutrinos will flow es- 
sentially in the radial direction, implying that the 
angular distribution becomes a delta- function in 
the forward direction so that (e)fl ux = (e). How- 
ever, in the trapping regions the two averages are 
very different because the distribution function is 
dominated by its isotropic term. 

To characterize the spectrum beyond the mean 
energy one can consider a series of moments (e n ) 
(Janka & Hillebrandt 1989a); we usually limit our- 
selves to n = 1 and 2. Note that a Fermi-Dirac 
distribution at zero chemical potential yields 



if!) 

(e}2 



486000 CsCs 

49 7T 8 



1.3029. 



(8) 



For a Maxwell-Boltzmann distribution this quan- 
tity would be 4/3. Following Raffelt (2001) we 
further define the "pinching parameter" 



P - a{e)> 



(9) 



where p = 1 signifies that the spectrum is thermal 
up to its second moment, while p < 1 signifies 
a pinched spectrum (high-energy tail suppressed) , 
p > 1 an anti-pinched spectrum (high-energy tail 
enhanced) . An analogous definition applies to the 
pinching parameter ph ux of the flux spectrum by 
replacing (•) with (-)fl U x- 



8 



In some publications the root-mean-square en- 
ergy (e)rms is given instead of the average energy. 
The definition corresponding to Eq. (4) is 



(e)r 



/ OO &/^ 1 1 d/xe 3 /(e,/z) 
\ Io°° de I-i d ^ e f( e ^) 




(10) 

This characteristic spectral energy is useful for es- 
timating the energy transfer from neutrinos to the 
stellar medium in reactions with cross sections 
proportional to e 2 . For thermal neutrinos with 
vanishing chemical potential we find 



(^)rms 



930 
441 



7T T« 4.5622 T. 



(11) 



With Eq. (6) this corresponds to (e) w 0.691 (e) rms . 

Beyond the energy moments (e n ) and related 
parameters it is often useful to approximate the 
neutrino spectrum by a simple analytic fit. If one 
uses two parameters beyond the overall normal- 
ization one can adjust the fit to reproduce two 
moments, for example (e) and (e 2 ). In the lit- 
erature one frequently encounters an approxima- 
tion in terms of a nominal Fermi-Dirac distribu- 
tion characterized by a temperature T and a de- 
generacy parameter rj according to 



fv( e ) = 7- Tl \ 

1 + exp (|r - rj) 



(12) 



(Janka & Hillcbrandt 1989a). In Fig. 7 we show 
(e/T) and p as a function of 77. Up to second order, 
expansions are 

(e/T) w 3.1514 + 0.1250 77 + 0.0429 rf 

p w 1 -0.0174 77- 0.0046 rj 2 . (13) 

These expansions are shown in Fig. 7 as dashed 
lines. 

Using a nominal Fermi-Dirac fit to approximate 
the spectrum is physically motivated because a 
truly thermal neutrino flux would follow this be- 
havior. On the other hand, the neutrino flux emit- 
ted from a SN core is not very close to being ther- 
mal so that the limiting behavior of the fit function 
is not a strong argument. Therefore, we consider 
an alternative fit function for which analytic sim- 
plicity is the main motivation, 



Me) = (§)". 



-(a+l)e/e 



(14) 



A 
H 

u 
V 
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Fig. 7. — Mean energy and pinching parameter 
as a function of the degeneracy parameter for a 
Fermi-Dirac distribution. As dashed lines we show 
the expansions given in Eq. (13). 

For any value of a we have (e) = e while 



if!) 

(e) 2 



2 + a 
T+~a 



(15) 



Put another way, e is the average energy while a 
represents the amount of spectral pinching. For 
general moments the analogous relation is 



(e fc ) 
{e k-i } 



a 



1 



(16) 



In the upper panel of Fig. 8 we show f a (e) , the 
integral normalized to unity, for several values of 
a. The broadest curve is for a — 2 while for the 
narrower ones the width 



w 



= V(t 2 ) - <c> 2 



(17) 



was decreased in 10% decrements as shown in Ta- 
ble 2. The middle panel of Fig. 8 shows the cor- 
responding curves f v (e) with the rj- values given in 
Table 2. The broadest curves in each panel are 
identical and correspond to e 2 exp(— 3e/e) with a 
width wq — (e) / \/3. 
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Table 2 

Parameters for fit-functions of Fig. 8. 



Width 


a 


V 


wo = (e)/\/3 


2. 


— oo 


0.9 w 


2.7037 


1.1340 


0.8 w 


3.6875 


2.7054 


0.7w a 


5.1225 


4.4014 


0.6 u>o 


7.3333 


6.9691 


0.5 w 


11. 


13.892 



The limiting behavior of j a {e) for large a is 
S(e — e) while for f v (e) the limiting width is 
Wo/y/5 ~ 0.44721 wo- Evidently the curves / a (e) 
can accommodate a much broader range of widths 
than the curves /^(e). 

We will find that the neutrino spectra are al- 
ways fit with parameters in the range 2 $ a < 5 or 
< 77 $ 4, i.e. with a width above about 0.75 wq- 
In the bottom panel of Fig. 8 we show the ratios 
of the fit functions for the widths down to 0.7 wo- 
Except for the lowest energies and very high ener- 
gies the two fit functions are equivalent to better 
than 10%. Therefore, the two types of fits are 
largely equivalent for most practical purposes. 

On the basis of a few high-statistics Monte- 
Carlo runs we will show in Sec. 3.5 that the numer- 
ical spectra are actually better approximated over 
a broader range of energies by the "power-law" fit 
functions / Q (e). In addition, these functions are 
more flexible at representing the high-energy tail 
of the spectrum that is most relevant for studying 
the Earth effect in neutrino oscillations. 

3.2. Monte Carlo Set Up 

We have run our Monte Carlo code, that is 
described in Appendix A, for the stellar back- 
ground models introduced in Sec. 2.2 and for dif- 
ferent combinations of energy-changing neutrino 
reactions. Our main interest is to assess the im- 
pact of the scattering atmosphere on the flux and 
spectrum formation. Therefore, it is sufficient to 
simulate the neutrino transport above some radius 
where we have to specify a boundary condition. 

We always use a blackbody boundary condition 
at the bottom of the atmosphere, i.e. we assume 



neutrinos to be in LTE at the local temperature 
and the appropriate chemical potential; for and 
v T the latter is taken to vanish. As a consequence 
of this boundary condition, the luminosity emerg- 
ing at the surface is generated within the computa- 
tional domain and calculated by our Monte Carlo 
transport. A small flux across the inner boundary 
develops because of the negative gradients of tem- 
perature and density in the atmosphere, but its 
magnitude depends on the radial resolution of the 
neutron-star atmosphere and will not in general 
correspond to the physical diffusive flux. But as 
long as the flux is small compared to the luminos- 
ity at the surface, the emerging neutrino spectra 
will not depend on the lower boundary condition. 
Usually it is sufficient to place the inner grid radius 
deeper in the star than the thermalization depth 
of the dominant pair process. 

The shallow energy dependence of the thermal- 
ization depth of the nucleon bremsstrahlung im- 
plies that whenever we include this process it is 
not difficult to choose a reasonable location for the 
lower boundary. Taking the latter too deep in the 
star is very CPU-expensive as one spends most of 
the simulation for calculating frequent scatterings 
of neutrinos that are essentially in LTE. 

We always include scattering as the main 
opacity source. For energy exchange, we switch on 
or off bremsstrahlung (b), nucleon recoil (r), scat- 
tering on electrons (s), e+e~ pair annihilation (p), 
and v c D c annihilation (n). We never include inelas- 
tic nucleon scattering v M NN — > NNf M as this pro- 
cess is never important relative to recoil (Raffclt 
2001). Likewise, we ignore scattering on v e and v c 
which is always unimportant if f M e± i s included 
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12 3 



e/{e) 




Q L. I I I I I I I I I I I I I I 

1 2 3 

e/{e) 

Fig. 8. — Normalized fit functions. Upper panel: 
"Power law" according to Eq. (14). Middle panel: 
Fermi-Dirac fit according to Eq. (12). In both 
panels the broadest curve corresponds to /(e) = 
e 2 exp(— 3e/e), i.e. to a = 2 and rj = — oo, re- 
spectively. For the other curves the width was 
decreased in decrements of 10%, see Table 2. Bot- 
tom panel: Ratio of the fits f a /f v for the first four 
cases. 

(Buras et al. 2002). We also neglect v^v^ or v^v^ 
scattering even though such processes may have 
a larger rate than some of the included leptonic 
processes. Processes of this type do not exchange 



energy between the neutrinos and the background 
medium. They are therefore not expected to affect 
the emerging fluxes and should also have a minor 
effect on the emitted spectra. 

3.3. Importance of Different Processes 

3.3.1. Accretion-Phase Model I 

Our first goal is to assess the relative impor- 
tance of different energy-exchange processes for 
the transport. As a first example we begin with 
our Accretion-Phase Model I. The results from our 
numerical runs are summarized in Table 3 where 
for each run we give (e)fl ux , our fit parameter 
a determined by Eq. (15), and the pinching pa- 
rameter pfl ux for the emerging flux spectrum, the 
temperature and degeneracy parameter of an ef- 
fective Fermi-Dirac spectrum producing the same 
first two energy moments, and the luminosity. 

The first row contains the muon neutrino flux 
characteristics of the original Boltzmann trans- 
port calculation by Messer. To make a connec- 
tion to these results we ran our code with the 
same input physics, i.e. i/ A ,e ± scattering (s) and 
e + e~ annihilation (p). There remain small differ- 
ences between the original spectral characteristics 
and ours. These can be caused by differences in 
the implementation of the neutrino processes, by 
the limited number of energy and angular bins in 
the Boltzmann solver, the coarser resolution of the 
radial grid in our Monte Carlo runs, and by our 
simple blackbody lower boundary condition. We 
interpret the first two rows of Table 3 as agreeing 
sufficiently well with each other that a detailed 
understanding of the differences is not warranted. 
Henceforth we will only discuss differential effects 
within our own implementation. 

In the next row (bsp) we include nucleon brems- 
strahlung which has the effect of increasing the 
luminosity by a sizable amount without affect- 
ing much the spectral shape. This suggests that 
bremsstrahlung is important as a source for v^v^ 
pairs, but that the spectrum is then shaped by the 
energy-exchange in scattering with e ± . In the next 
row we switch off e ± scattering (bp) so that no 
energy is exchanged except by pair-producing pro- 
cesses. The spectral energy indeed increases sig- 
nificantly However, the biggest energy-exchange 
effect in the scattering regime is nucleon recoil. 
In the next two rows we include recoil (brp) and 
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Table 3 

Monte Carlo results for Accretion-Phase Model I. 



Energy 


exc 


han 


ge 


(e)flux 


/ 2\ 

(e /flux 


a 


Pflux 


T 


V 


L v 


original 


run 




17.5 


388. 


2.7 


0.97 


5.2 


1.1 


14.4 




s 


P 




16.6 


362. 


2.2 


1.01 


5.3 


-0.3 


15.8 


b 


s 


P 




16.3 


351. 


2.1 


1.02 


5.4 


-2.2 


19.1 


b 




P 




17.8 


419. 


2.1 


1.02 


5.9 


-1.9 


20.1 


b r 




P 




15.1 


285. 


3.0 


0.96 


4.3 


1.6 


18.6 


b r 


s 


P 




14.2 


255. 


2.8 


0.98 


4.2 


1.1 


14.8 


b r 


s 


P 


n 


14.4 


264. 


2.7 


0.97 


4.3 


1.2 


17.6 


b 


s 


P 


n 


16.6 


358. 


2.3 


1.00 


5.2 


0.2 


21.7 




s 


P 


n 


16.9 


369. 


2.4 


0.99 


5.3 


0.4 


20.2 


b r 


s 






14.0 


251. 


2.6 


0.99 


4.3 


0.6 


13.1 


b r 


s 




n 


14.4 


263. 


2.7 


0.97 


4.3 


1.2 


17.0 


r 


s 


P 




14.5 


265. 


2.8 


0.97 


4.3 


1.2 


13.0 


r 


s 


P 


n 


14.7 


269. 


3.1 


0.96 


4.2 


1.7 


16.8 


b r 


sn 


P 


n 


14.3 


260. 


2.7 


0.97 


4.3 


1.2 


17.9 



Note. — For energy exchange, "b" refers to bremsstrahlung, "r" to 
recoil, "s" to scattering on electrons and positrons, "p" to e + e _ anni- 
hilation, "n" to v c v c annihilation, and "sn" to scattering on both e^ 
and v c , v c . We give (e)fl ux and T in MeV, (e 2 )fl ux in MeV 2 , and L v in 
10 51 erg s" 1 . 
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then additionally c ± scattering (brsp) , both lower- 
ing the spectral energies and also the luminosities. 

The picture of all relevant processes is com- 
pleted by adding v c v c pair annihilation (brspn), 
which is similar to e + e~ pair annihilation, but a 
factor of 2-3 more important (Buras et al. 2002). 
The luminosity is again increased, an effect which 
is understood in terms of our blackbody picture 
for the number and energy spheres. In the lower 
panel of Fig. 3 we see that i?therm moves to larger 
radii once "n" is switched on, the radiating surface 
of the "blackbody" increases and more pairs are 
emitted. For both "p" and "n" i?thcrm is strongly 
energy dependent and therefore it is impossible to 
define a sharp thermalization radius. 

Switching off "r" again (bspn) shows that also 
with "n" included, "r" really dominates the mean 
energy and shaping of the spectrum. 

To study the relative importance of the differ- 
ent pair processes, we switch off the leptonic ones 
(row "brs" ) and compare this to only the leptonic 
processes (row "rspn"). In this stellar model both 
types contribute significantly. Comparing then 
"brsp" with "brsn" shows that among the leptonic 
processes "n" is clearly more important than "p" . 

The last row "brsnpn" includes in addition to 
all other processes scattering on v e and v e . It was 
already shown by Buras et al. (2002) that this pro- 
cess is about half as important as scattering on e ± 
and its influence on the neutrino flux and spectra 
is negligible. We show this case for completeness 
but do not include scattering on v c and v e for any 
of our further models. 

In order to illustrate some of the cases of Ta- 
ble 3 we show in the upper panel of Fig. 9 several 
flux spectra from high-statistics Monte-Carlo runs. 
Starting again with the input physics of the origi- 
nal hydrodynamic simulation (sp) we add brems- 
strahlung (b), recoil (r), and finally v c v c pair an- 
nihilation (n). Each of these processes has a sig- 
nificant and clearly visible influence on the curves. 
The pair-creation processes ( "b" and "n" ) hardly 
change the spectral shape but increase the num- 
ber flux, whereas recoil (r) strongly modifies the 
spectral shape. In the lower panel of Fig. 9 we 
show the same curves, normalized to equal parti- 
cle fluxes. In this representation it is particularly 
obvious that the pair processes do not affect the 
spectral shape. 
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Fig. 9. — High-statistics spectra for Accretion- 
Phase Model I with different input physics as in 
Table 3. Upper panel: Differential particle fluxes. 
Lower panel: Spectra normalized to equal particle 
fluxes. 

The very different impact of pair processes and 
nucleon recoils has a simple explanation. The 
thermalization depth for the pair processes is 
deeper than that of the energy-exchanging re- 
actions, i.e. the "number sphere" is below the 
"energy sphere." Therefore, the particle flux is 
fixed more deeply in the star while the spectra are 
still modified by energy-exchanging reactions in 
the scattering atmosphere. 

3.3.2. Steep Power Law 

As another example we study the steep power- 
law model defined in Eq. (3) and Table 1. This 
model is supposed to represent the outer layers 
of a late-time proto-neutron star but without be- 
ing hydrostatically self-consistent. It connects di- 
rectly with Raffelt (2001), where the same pro- 
file was used in a plane parallel setup, studying 
bremsstrahlung and nucleon recoil. The results of 
our runs are displayed in Table 4 and agree very 
nicely with those obtained by Raffelt (2001), cor- 
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Table 4 

Monte Carlo results for the steep power-law model. 



Energy 


ex 


change 




Y c 


(e)flux 


(e 2 )fiux 


a 


Mux 


T 


r\ 


L v 













9^ S 
ZD.o 


G£9 


1 9 
1 .Z 


i.ii 






91 n 


D 


r 










A Q.7 


9 P. 
Z.D 


U.ao 


o.u 


u. / 


1 A K 
14.0 


D 




P 




u 


ZD. 4: 


son 


1.0 


1 .uo 






9^ R 


b 


I 


p 


_ 


0.05 


25.6 


908. 


1.6 


1.06 


Z 


Z 


23.2 


b 




P 




0.5 


25.5 


917. 


1.4 


1.08 






21.6 


b 




s p 







24.2 


787. 


1.9 


1.03 






24.5 


b 




s p 




0.05 


23.8 


753. 


2.0 


1.02 






24.5 


b 




s p 




0.5 


21.3 


591. 


2.3 


1.00 


6.8 


-0.3 


23.1 


b 


r 


P 




0.5 


20.0 


507. 


2.7 


0.98 


6.0 


1.0 


16.8 


b 


r 


s p 







20.3 


518. 


2.9 


0.97 


5.9 


1.4 


19.7 


b 


r 


s p 




0.05 


20.3 


518. 


2.9 


0.97 


5.9 


1.4 


19.5 


b 


r 


s p 




0.5 


19.6 


488. 


2.7 


0.98 


5.9 


1.1 


18.7 


b 


r 


s p 


n 





20.7 


535. 


3.0 


0.96 


5.8 


1.8 


23.9 


bx3 


r 


s p 


n 


0.05 


20.3 


522. 


2.7 


0.97 


6.0 


1.3 


24.2 


b 


r 


s p 


n 


0.05 


20.6 


530. 


3.0 


0.96 


5.9 


1.7 


23.8 


bx0.3 


r 


s p 


n 


0.05 


20.7 


534. 


3.1 


0.96 


5.8 


1.8 


23.4 


b 


r 


s p 


n 


0.5 


19.8 


499. 


2.7 


0.97 


5.9 


1.2 


21.4 



Table 5 

Monte Carlo results for the shallow power-law model. 



Energy exchange 


Y c 


(e)flux 


(e 2 )flux 


a 


Mux 


T 


V 


L v 


b 










27.7 


1120. 


1.2 


1.12 






20.3 


b r 










20.1 


521. 


2.5 


0.99 


6.3 


0.4 


13.4 


b 




p 







27.7 


974. 


2.7 


0.98 


8.3 


1.0 


43.1 


b 




p 




0.05 


27.9 


990. 


2.7 


0.98 


8.3 


1.1 


43.3 


b 




p 




0.5 


28.3 


1019. 


2.7 


0.98 


8.5 


1.0 


38.3 


b 


s 


p 







25.5 


830. 


2.6 


0.98 


7.6 


1.1 


46.2 


b 


s 


p 




0.05 


25.4 


815. 


2.8 


0.97 


7.5 


1.2 


46.3 


b 


s 


p 




0.5 


23.5 


706. 


2.6 


0.98 


7.1 


1.0 


44.8 


b r 




p 




0.5 


22.5 


624. 


3.3 


0.95 


6.1 


2.2 


33.1 


b r 


s 


p 







22.3 


612. 


3.3 


0.95 


6.1 


2.1 


39.6 


b r 


s 


p 




0.05 


22.2 


609. 


3.2 


0.95 


6.1 


2.1 


39.1 


b r 


s 


p 




0.5 


21.7 


585. 


3.1 


0.95 


6.1 


1.9 


39.2 


b r 


s 


p 


n 





22.2 


608. 


3.3 


0.94 


6.0 


2.2 


54.7 


b r 


s 


p 


n 


0.05 


22.4 


615. 


3.4 


0.94 


6.1 


2.3 


54.9 


b r 


s 


p 


n 


0.5 


21.8 


587. 


3.3 


0.95 


6.1 


1.9 


51.3 
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responding to our cases "b" and "br" . 

For investigating the importance of leptonic 
processes, we run our code with a variety of neu- 
trino interactions and in addition assume a con- 
stant electron fraction Y e throughout the whole 
stellar atmosphere. This assumption is somewhat 
artificial, but gives us the opportunity to study 
extreme cases in a controlled way. In the relevant 
region Y c = 0.5 yields the highest possible electron 
density. In addition we study the electron fraction 
being one order of magnitude smaller, Y c = 0.05, 
and finally the extreme case with an equal number 
of electrons and positrons, Y c = 0. 

The first leptonic process we consider is c + e~ 
pair annihilation. Comparing the rows "bp" with 
the row "b" shows a negligible effect on the spec- 
trum, but a rise in luminosity. Increasing Y e 
brings the luminosity almost back to the "b" 
case, because the electron degeneracy rises and the 
positron density decreases so that the pair process 
becomes less important. 

Adding scattering on e ± forces the transported 
neutrinos to stay closer to the medium tempera- 
ture, i.e. reduces their mean energy. Of course, the 
scattering rate increases with the number of elec- 
trons and positrons, i.e. for higher Y e we get lower 
spectral energies. For the luminosity the situation 
is more complicated. Since the neutrino flux en- 
ergies decrease when we switch on e^ scattering 
we would expect a lower luminosity. However, the 
opacity of the medium to neutrinos is strongly en- 
ergy dependent and low energy neutrinos can es- 
cape more easily than high-energy ones, increasing 
the number flux. On balance, the "bsp" luminosi- 
ties are larger compared to the "bp" ones. 

To compare the scattering on e 1 * 1 with that on 
nucleons, we turn off "s" again and instead switch 
on recoil (r) . Qualitatively, the energy exchange is 
very different from the earlier case. In the scatter- 
ing on a neutrino can exchange a large amount 
of energy, while for scattering on nucleons the en- 
ergy exchange is small. But since neutrino-nucleon 
scattering is the dominant source of opacity that 
keeps the neutrinos inside the star, the scatter- 
ings are very frequent. This leads to a stronger 
suppression in the high-energy tail of the neutrino 
spectrum and therefore to a visibly smaller mean 
flux energy and lower effective spectral tempera- 
ture, but higher effective degeneracy. Many nu- 
cleoli scatterings, however, are needed to down- 



grade the high-energy neutrinos (different from 
e"* 1 scattering). Therefore neutrinos stay longer at 
high energies and experience a larger opacity and a 
larger amount of backscattering. This suppresses 
the neutrino flux significantly. 

In the runs including both scattering reactions 
(brsp), we find a mixture of the effects of e ± and 
nucleon scatterings and an enhanced reduction of 
the mean flux energy. 

Finally, adding the neutrino pair process yields 
almost no change in energy and pinching, but an 
increased luminosity as expected from the analo- 
gous case in Sec. 3.3.1. Although this profile is 
rather steep, leptonic pair processes are still im- 
portant (Fig. 5). 

In order to estimate the sensitivity to the exact 
treatment of nucleon brcmsstrahlung we have per- 
formed one run with the bremsstrahlung rate arti- 
ficially enhanced by a factor of 3, and one where it 
was decreased by a factor 0.3. All other processes 
were included. The emerging fluxes and spec- 
tra indeed do not depend sensitively on the exact 
strength of bremsstrahlung as argued in Sec. 2.1. 

3.3.3. Shallow Power Law 

For the shallow power law almost the same dis- 
cussion as for the steep case applies. As we can 
already infer from Fig. 6, leptonic processes are 
more important. This leads to a much higher in- 
crease of the neutrino flux once "p" or "n" are 
included, and to stronger spectral pinching when 
c+e~ annihilation is switched on. Scattering on 
c ± downgrades the transported neutrino flux by a 
larger amount. 

3.4. The Effect of Binning 

Evidently nucleon recoil plays an important role 
for the l/fj, spectrum formation. It is straight- 
forward to implement this reaction in our Monte 
Carlo approach, but it may be more difficult in 
those treatments of neutrino transport that rely on 
binned energy spectra. The energy exchange in a 
given i/N collision is relatively small (see e.g. Raf- 
felt 2001) so that one may lose this effect if the 
spectrum is too crudely binned. To test the im- 
pact of binning we have performed two runs for the 
Accretion Phase Model I where we fix the neutrino 
energies to a small number of values. After every 
interaction the final-state energy is set to the cen- 
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tral value of the energy bin it falls into. We have 
chosen the 17 logarithmically spaced bins on the 
interval from to around 380 MeV that are used 
in the simulations of the Garching SN group. 

In Fig. 10 we compare the emergent flux spectra 
from runs with binning (histograms) and without 
(smooth curves). For both cases recoil was once 
included (brspn) and once not (bspn). The mean 
energies calculated from either binned or unbinned 
runs agree to better than 0.5%, and also the shapes 
are well reproduced, although there are slight dif- 
ferences around the spectral peak. We conclude 
that the impact of recoil is well accounted for in 
our runs with discrete energies. Therefore, the en- 
ergy grid of the Garching group should well suffice 
to reproduce the main spectral impact of nucleon 
recoil. 
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Fig. 10. — Comparison of runs with binned and 
unbinned neutrino energies. 



3.5. Detailed Spectral Shape 

Thus far we have characterized the neutrino 
spectra by a few simple parameters. However, it 
is extremely useful to have a simple analytic fit to 
the overall spectrum that can be used, for exam- 
ple, to simulate the response of a neutrino detector 
to a SN signal. To study the quality of different fit 
functions we have performed a few high-statistics 
Monte-Carlo runs for the Accretion Phase Model I, 
including all interaction processes. Moreover, we 
have performed these runs for the flavors v Cl D e , 
and v^. (A detailed discussion of the flavor de- 
pendence is deferred to Sec. 4.) 

In order to get smooth spectral curves we have 
averaged the output of 70,000 time steps. In addi- 



tion we have refined the energy grid of the neutrino 
interaction rates. Both measures leave the previ- 
ous results unaffected but increase computing time 
and demand for memory significantly. 

In Fig. 11 we show our high-statistics Monte 
Carlo (MC) spectra together with the a-fit func- 
tion f a (e) defined in Eq. (14) and the ?7-fit func- 
tion f v (c) of Eq. (12). The analytic functions can 
only fit the spectrum well over a certain range of 
energies. We have chosen to optimize the fit for 
the event spectrum in a detector, assuming the 
cross section scales with e 2 . Therefore, we actu- 
ally show the neutrino flux spectra multiplied with 
e 2 . Accordingly, the parameters a and e, as well 
as r\ and T and the normalizations are determined 
such that the energy moments (e 2 ), (e 3 ), and (e 4 ) 
are reproduced by the fits. 

Below each spectrum we show the ratio of our 
MC results with the fit functions. In the energy 
range where the statistics in a detector would be 
reasonable for a galactic SN, say from 5-10 MeV 
up to around 40 MeV, both types of fits represent 
the MC results nicely. However, in all cases the 
a-fit works somewhat better than the 7/-fit. 

We have repeated this exercise for the steep 
power-law models with q = 2.5 and the one with 
q = 3.5. The quality of the fits is comparable to 
the previous example. 

3.6. Summary 

We find that the spectra are reasonably well 
described by the simple picture of a blackbody 
sphere determined by the thcrmalization depth of 
the nucleonic bremsstrahlung process, the "filter 
effect" of the scattering atmosphere, and energy 
transfers by nucleon recoils. This is also true for 
the flux in case of steep neutron-star atmo- 
spheres. For more shallow atmospheres pair anni- 
hilation (e + e _ and v e P e ), however, yields a large 
contribution to the emitted flux and e scatter- 
ing reduces the mean flux energy significantly. It is 
therefore important for state-of-the art transport 
calculations to include these leptonic processes. 
The traditional process e + e~ — > v^v^ is subdomi- 
nant compared to v c v c — > v^v^ as previously found 
by Buras et al. (2002). The relative importance of 
the various reactions depends on the stellar profile. 

Neutrinos emitted from a blackbody surface 
and filtered by a scattering atmosphere with- 
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out recoils and leptonic processes have an anti- 
pinched spectrum (Raffelt 2001). However, af- 
ter all energy-exchanging reactions have been in- 
cluded we find that the spectra are always pinched. 
When described by effective Fermi-Dirac distribu- 
tions, the nominal degeneracy parameter r\ is typ- 
ically in the range 1-2, depending on the profile 
and electron concentration. 
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Fig. 11. — High-statistics spectra for Accretion-Phase Model I including all interaction processes. The Monte 
Carlo (MC) results are shown as crosses, the analytic fit functions as smooth lines. The left-hand panels use 
as fits / a (e) according to Eq. (14), the right-hand panels f v (e) according to Eq. (12). Below the spectra we 
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Table 6 

Comparing Monte Carlo results for different flavors. 



Model and Flavor 




(e)flux 


(e 2 )flux 


\£)fiux 


a 


Pflux 


T 


V 




Accretion-Phase Model I 


















Original 
























17.5 


388. 


1.20 


2.7 


0.97 


5.2 


1.1 


14.4 






14.6 


253. 


1 


4.4 


0.91 


3.5 


3.4 


29.2 






12.5 


190. 


0.86 


3.6 


0.93 


3.2 


2.8 


30.8 


Our runs 




















*V p m (" S P") 




16.6 


362. 


1.19 


2.2 


1.01 


5.3 


-0.3 


15.8 


P M ("brspn") 




14.3 


260. 


1.02 


2.7 


0.97 


4.3 


1.2 


17.9 


(weak magnetism) 




14.2 


254. 


1.01 


2.9 


0.96 


4.1 


1.6 


17.4 


£^ (weak magnetism) 




14.9 


281. 


1.06 


2.8 


0.97 


4.3 


1.5 


18.3 






14.0 


237. 


1 


3.8 


0.93 


3.6 


2.7 


31.7 






11.8 


175. 


0.84 


2.9 


0.97 


3.4 


1.4 


31.9 


Accretion-Phase Model II 


















Original 




















^, ^ 




17.2 


380. 


1.09 


2.5 


0.98 


5.2 


0.8 


32.4 






15.8 


300. 


1 


4.0 


0.92 


4.0 


3.0 


68.1 






12.9 


207. 


0.82 


3.1 


0.96 


3.7 


1.7 


65.6 


Our runs 
























15.7 


317. 


1.02 


2.5 


0.98 


4.8 


0.8 


27.8 






15.4 


283. 


1 


4.2 


0.92 


3.8 


3.2 


73.5 






13.0 


207. 


0.84 


3.4 


0.95 


3.6 


2.1 


73.9 


Steep Power Law p — 10 


















q = 2.5 




















^, ^ 


0.15 


20.4 


525. 


1.10 


2.8 


0.96 


5.9 


1.5 


23.5 


v c 


0.15 


18.5 


413. 


1 


3.8 


0.92 


4.6 


3.0 


23.5 




0.15 


12.7 


198. 


0.69 


3.4 


0.94 


3.4 


2.4 


12.8 




0.2 


20.4 


521. 


1.14 


3.0 


0.97 


5.9 


1.5 


23.3 




0.2 


17.9 


383. 


1 


4.1 


0.92 


4.4 


3.1 


11.7 


v c 


0.2 


13.4 


218. 


0.75 


3.7 


0.93 


3.4 


2.9 


24.4 


q = 3.0 






















0.1 


17.7 


393. 


1.14 


2.9 


0.96 


5.0 


1.8 


12.7 




0.1 


15.5 


289. 


1 


3.9 


0.93 


4.0 


2.8 


8.8 


v c 


0.1 


10.5 


132. 


0.68 


4.1 


0.92 


2.6 


3.0 


6.6 


q = 3.5 




















^, ^ 


0.07 


15.8 


310. 


1.22 


3.1 


0.95 


4.4 


2.1 


7.9 


^ M (weak magnetism) 


0.07 


15.5 


296. 


1.19 


3.3 


0.94 


4.2 


2.3 


7.7 


(weak magnetism) 


0.07 


16.5 


337. 


1.27 


3.2 


0.95 


4.5 


2.1 


8.3 


v c 


0.07 


13.0 


207. 


1 


3.4 


0.94 


3.5 


2.3 


4.3 


v a 


0.07 


9.4 


103. 


0.72 


5.0 


0.90 


2.1 


3.9 


4.1 


Shallow Power Law p — 


5, q — 


1 


















0.3 


22.0 


596. 


1.14 


3.3 


0.94 


6.0 


2.2 


53.9 




0.3 


19.3 


440. 


1 


4.5 


0.91 


4.5 


3.7 


85.7 



Table 6 — Continued 



Model and Flavor Y c (e)fl ux 


(e 2 )fl ux 


(e)flux 
7 — : a p flux 


T 


V L v 










u e 0.3 14.7 


262. 


0.76 3.7 0.93 


3.8 


2.7 56.5 



Note. — We give (e) flux and T in MeV, (e 2 ) flux in MeV 2 , L v in 10 51 erg s -1 . 
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4. COMPARING DIFFERENT 
FLAVORS 

4.1. Monte Carlo Study 

The new energy-exchange channels studied in 
the previous section lower the average energies. 
In order to compare the fluxes and spectra with 
those of v c and D c we perform a new series of runs 
where we include the full set of relevant micro- 
physics for Vy, and also simulate the transport of 
v e and v c . 

The microphysics for the interactions of v e and 
D c is the same as in Janka & Hillebrandt (1989a, b), 
i.e. charged-current reactions of with nucleons, 
iso-energetic scattering on nucleons, scattering on 
e ± , and e + e~ pair annihilation. In principle one 
should also include nucleon bremsstrahlung and 
the effect of nucleon recoils for the transport of v e 
and D e , but their effects will be minimal. There- 
fore, we preferred to leave the original working 
code unmodified for these flavors. 

In the first three rows of Table 6 we give the 
spectral characteristics for the Accretion-Phase 
Model I from the original simulation of Messer. 
The usual hierarchy of average neutrino energies 
is found, i.e. (e v .) : (e 5 .) : (e V(t ) = 0.86 : 1 : 1.20. 
The luminosities are essentially equal between v e 
and v e while v^, D^, v r , and v T each provide about 
half of the v c luminosity. 

Our Monte Carlo runs of this profile establish 
the same picture for the same input physics. Al- 
though our mean energies are slightly offset to 
lower values for all flavors relative to the origi- 
nal run, our energies relative to each other are 
(e„.) : (e 5 .) : (e„„) = 0.84 : 1 : 1.19 and thus very 
similar. However, once we include all energy ex- 
changing processes we find 0.84 : 1 : 1.02 instead. 
Therefore, (e„ ) no longer exceeds (e Pe ) by much. 
The luminosity of is about half that of v c or 
D e which are approximately equal, in rough agree- 
ment with the original results. Even though the 
additional processes lower the mean energy of 
they yield a more than 10% higher luminosity, 
mainly due to v e v e annihilation. 

As another example of an accreting proto- 
neutron star we use the Accretion-Phase Model II. 
The neutrino interactions included in this model 
were nucleon bremsstrahlung, scattering on c^, 
and e + e~ annihilation. Nucleon correlations, ef- 



fective mass, and recoil were taken into account, 
following Burrows & Sawyer (1998, 1999), as 
well as weak magnetism effects (Horowitz 2002) 
and quenching of qa at high densities (Carter & 
Prakash 2002). All these improvements to the 
traditional microphysics affect mainly v^, and to 
some degree also v c . Weak magnetism terms de- 
crease the nucleon scattering cross sections for 
more strongly than they modify scatterings. In 
this hydrodynamic calculation, however, and 
Vy were treated identically by using the average 
of the corresponding reaction cross sections. The 
effects of weak magnetism on the transport of 
and Vy are therefore not included to very high 
accuracy. Note, moreover, that the original data 
come from a general relativistic hydrodynamic 
simulation with the solution of the Boltzmann 
equation for neutrino transport calculated in the 
comoving frame of the stellar fluid. Therefore the 
neutrino results are affected by gravitational red- 
shift and, depending on where they are measured, 
may also be blucshiftcd by Doppler effects due to 
the accretion flow to the nascent neutron star. 

Our Monte Carlo simulation in contrast was 
performed on a static background without general 
relativistic corrections. It includes bremsstrah- 
lung, recoil, e + e~ pair annihilation, scattering on 

, and v c v c annihilation, i.e. our microphysics is 
similar but not identical with that used in the orig- 
inal run. As an outer radius we took 100 km; all 
flux parameters are measured at this radius be- 
cause farther out Doppler effects of the original 
model would make it difficult to compare the re- 
sults. Keeping in mind that we use very differ- 
ent numerical approaches and somewhat different 
input physics, the agreement in particular for v c 
and v c is remarkably good. This agreement shows 
once more that our Monte Carlo approach likely 
captures at least the differential effects of the new 
microphysics in a satisfactory manner. 

In Fig. 12 we compare our calculations for the 
Accretion-Phase Model II with those of the orig- 
inal simulation. The step-like curve again repre- 
sents the mean energy of neutrinos in LTE for zero 
chemical potential. The smooth solid line is the 
mean energy (e) from our runs, the dotted (lower) 
line gives (e)fl ux . The crosses are the correspond- 
ing results from the original runs. For the trans- 
port of Vy our inner boundary is R m = 16 km, 
while for v c and v e we use R m = 24 km. For v c 
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and P e the charged-current processes (urea) keep 
these neutrinos in LTE up to larger radii than pair 
processes in the case of f M . With our choice of i?; n 
the neutrinos are in LTE within the innermost ra- 
dial zones. 

The results are similar to the Accretion-Phase 
Model I. The luminosities are not equipartitioned 
but instead follow roughly L Vo w L 5e w 2L Ufi . 
The ratios of mean energies arc {e u J : (e Po ) : 
(e„ ) — 0.82 : 1 : 1.09 in the original run and 
0.84 : 1 : 1.02 in our run. 

In summary, both accretion-phase models agree 
reasonably well in the (e„ e ) : (e Po ) ratio for all runs. 
Moreover, using traditional input physics one finds 
something like (e Po ) : (e„ ) = 1 : 1.20. Depending 
on the implementation of the new input physics 
and depending on the model one finds results be- 
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Fig. 12. — Comparison of the Accretion-Phase 
Model II calculations. Continuous lines show our 
Monte Carlo runs while crosses represent the orig- 
inal simulation by Rampp. The steps correspond 
to (e) = 3.15 T. 



tween (e Pc ) : (e^) = 1 : 1.02 and 1 : 1.09. The 
higher ratio in Rampp's simulation could be due 
to the inclusion of weak magnetism which tends 
to raise (e„ ) more than (e Pc ). 

In order to estimate the corresponding results 
for later stages of the proto-neutron star evolution 
we employ our steep power-law model. We vary 
the power q of the temperature profile within a 
reasonable range so that q/p = 0.25-0.35, with q 
and p defined in Eq. (3). Y e is fixed by demanding 
roughly equal number fluxes for v e and D e because 
a few seconds after bounce deleptonization should 
be essentially complete. The fluxes of these neu- 
trinos depend very sensitively on Y c so that this 
constraint is only reached to within about 30% 
without tuning Y c to three decimal places. How- 
ever, the mean energies are rather insensitive to 
the exact value of Y c . This is illustrated by the 
steep power-law model with q = 2.5 where we 
show results for Y c = 0.15 and 0.20. The num- 
ber fluxes of u e and v c differ by less than 30% for 
Y c = 0.15, but differ by a factor of 3 for Y c = 0.2. 
At the same time, the average spectral energies 
barely change. 

The ratios of mean energies are not very differ- 
ent from those of the accretion-phase models. Of 
course, the absolute flux energies have no physical 
meaning because we adjusted the stellar profile in 
order to obtain realistic values. For the luminosi- 
ties we find L„ o < L Vli , different from the accretion 
phase. The steep power law implies that the ra- 
diating surfaces are similar for all flavors so that 
it is not surprising that the flavor with the largest 
energies also produces the largest luminosity. 

We find that (e v ) always exceeds (e Po ) by a 
small amount, the exact value depending on the 
stellar model. During the accretion phase the en- 
ergies seem to be almost identical, later they may 
differ by up to 20%. We have not found a model 
where the energies differ by the large amounts 
which are sometimes assumed in the literature. At 
late times when Y c is small the microphysics gov- 
erning D e transport is closer to that for than 
at early times. Therefore, one expects that at late 
times the behavior of v c is more similar to than 
at early times. We do not see any argument for 
expecting an extreme hierarchy of energies at late 
times for self-consistent stellar models. 

We never find exact equipartition of the flavor- 
dependent luminosities. Depending on the stellar 
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profile the fluxes can mutually differ by up to a 
factor of 2 in either direction. 

4.2. Weak Magnetism 

Weak magnetism causes a significant correction 
to the neutrino-nucleon cross section that arises 
due to the large anomalous magnetic moments 
of protons and neutrons (Vogel & Bcacom 1999, 
Horowitz 2002). It increases the neutrino interac- 
tion rate but lowers the rate for anti-neutrinos. It 
is expected to be a small correction in the SN con- 
text, but has never been implemented so far. Fol- 
lowing Horowitz (2002) we add weak magnetism 
to our nucleon-recoil rate as given in Appendix B. 

Our Monte Carlo code transports only one 
species of neutrinos at a time. In order to test 
the impact of weak magnetism we assumed that 
a chemical potential for would build up, and 
assumed a fixed value for the degeneracy pa- 
rameter throughout our stellar model. We then 
iterated several runs for and with different 
degeneracy parameters until their particle fluxes 
were equal because in a stationary state there will 
be no net flux of /i-lepton number. 

We performed this procedure for our Accretion- 
Phase Model I and our steepest power law; the re- 
sults are summarized in Table 6. In both cases the 
mean energies of go down by 2% and go up for 

by 4%. The mean luminosities arc unaffected. 
We conclude that weak-magnetism corrections are 
small. Transporting and separately in a 
self-consistent hydrodynamic simulation is prob- 
ably not worth the cost in computer time. 

4.3. Previous Literature 

There is a large recent body of literature quoted 
in our introduction where the effect of flavor oscil- 
lations on SN neutrino spectra and fluxes is stud- 
ied. Many of these papers assumed that (e„ ) is 
much larger than (ep e ) and that the luminosities 
between all flavors were exactly equipartitioned. 
Our findings here are almost orthogonal to this 
perception. Where does it come from? 

To the best of our knowledge, the microphysics 
employed for transport is roughly the same in 
all published simulations. It includes iso-energetic 
scattering on nuclcons, e + c~ annihilation and 
v^e^- scattering. Of course, the transport method 
and the numerical implementation of the neutrino 



processes differ in the codes of different groups. 
The new reactions and nucleon recoil lower (e^) 
and modify the luminosities, but not by such a 
large amount as to explain a completely different 
paradigm. Therefore, we have inspected the previ- 
ous literature and collect a representative sample 
of pertinent results in Table 7. Note that the sim- 
ulations discussed below did not in all cases use 
the same stellar models and equations of state for 
the dense matter in the supernova core. 

We begin with the simulations of the Livermore 
group who find robust explosions by virtue of the 
neutron-finger convection phenomenon. Neutrino 
transport is treated in the hydrodynamic models 
with a multigroup flux-limited diffusion scheme. 
Maylc, Wilson, & Schramm (1987) gave detailed 
results for their SN simulation of a 25 Mq star. 
For half a second after bounce they obtained a 
somewhat oscillatory behavior of the neutrino lu- 
minosities. After the prompt peak of the electron 
neutrino luminosity, they got L„ c w L Po w 2 L V)i w 
50-130 x 10 51 erg s _1 . After about one second the 
values stabilize. This calculation did not produce 
the "standard" hierarchy of energies. However, 
there is clearly a tendency that v c behave more 
similar to at late times. 

The most recent published Livermore simula- 
tion is a 20 Mq star (Totani et al. 1998). It shows 
an astonishing degree of luminosity equipartition 
from the accretion phase throughout the early 
Kelvin-Hclmholtz cooling phase. About two sec- 
onds after bounce the flux falls off more slowly 
than the other flavors. In Table 7 we show repre- 
sentative results for an early and a late time. The 
mean energies and their ratios are consistent with 
what we would have expected on the basis of our 
study. 

With a different numerical code, Bruenn (1987) 
found for a 25 Mq progenitor qualitatively differ- 
ent results for luminosities and energies. At about 
0.5 s after bounce the luminosities and energies 
became stable at the values given in Table 7. This 
simulation is an example for an extreme hierarchy 
of mean energies. 

In Burrows (1988) all luminosities are said to 
be equal. In addition it is stated that for the first 
5 seconds {t v ) w 24 MeV and the relation to the 
other flavors is (e„ e ) : (e Pc ) : (e^) = 0.9 : 1 : 1.8. 
Detailed results are only given for P e , so we are not 
able to add this reference to our table. The large 



23 



Table 7 

Flavor dependent flux characteristics from the literature. 





tpb 


<*,.> 


(<*.> 








) 








Mayle et al. (1987) 


1.0 


12 


24 


22 


0.50 


: 1 


: 0.92 


20 


20 


20 


Totani et al. (1998) 


0.3 


12 


15 


19 


0.80 


: 1 


: 1.26 


20 


20 


20 




10 


11 


20 


25 


0.55 


: 1 


: 1.25 


0.5 


0.5 


1 


Brucnn (1987) 


0.5 


10 


12 


25 


0.83 


: 1 


: 2.08 


3 


5 


16 


Myra & Burrows (1990) 


0.13 


11 


13 


24 


0.85 


: 1 


: 1.85 


30 


30 


16 


Janka & Hillebrandt (1989b) 


0.3 


8 


14 


16 


0.57 


: 1 


: 1.14 


30 


220 


65 


Suzuki (1990) 


1 


9.5 


13 


15 


0.73 


: 1 


: 1.15 


4 


4 


3 




20 


8 


10 


9 


0.80 


: 1 


: 0.90 


0.3 


0.3 


0.07 


Suzuki (1991) 


1 


9.5 


13 


15 


0.73 


: 1 


: 1.15 


3 


3 


3 




15 


8 


9 


9.5 


0.89 


: 1 


: 1.06 


0.4 


0.4 


0.3 


Suzuki (1993) 


1 


9 


12 


13 


0.75 


: 1 


: 1.08 


3 


3 


3 




15 


7 


8 


8 


0.88 


: 1 


: 1.00 


0.3 


0.3 


0.3 


Accretion-Phase Model I (original) 


0.32 


13 


15 


18 


0.86 


: 1 


: 1.20 


31 


29 


14 


Accretion-Phase Model I (our run) 


0.32 


12 


14 


14 


0.84 


: 1 


: 1.02 


32 


32 


18 


Accretion-Phase Model II (original) 


0.15 


13 


16 


17 


0.82 


: 1 


: 1.09 


66 


68 


32 


Accretion-Phase Model II (our run) 


0.15 


13 


15 


16 


0.84 


: 1 


: 1.02 


74 


74 


28 


Buras et al. (personal comm.) 


0.25 


14.1 


16.5 


16.8 


0.85 


: 1 


: 1.02 


43 


44 


32 


The following lines show (e) rms 


instead of (e) 


















Mezzacappa et al. (2001) 


0.5 


16 


19 


24 


0.84 


: 1 


: 1.26 


25 


25 


8 


Liebendorfer et al. (2001) 


0.5 


19 


21 


24 


0.90 


: 1 


: 1.14 


30 


30 


10 



Note. — We give the time post bounce (tpb) in s, (e) in MeV, and L v in 10 51 erg s 1 . 
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variety of models investigated by Burrows (1988) 
and the detailed results for i> e go beyond the scope 
of our brief description. In a later paper Myra & 
Burrows (1990) studied a 13 M Q progenitor model 
and found the extreme hierarchy of energies shown 
in our table. 

With the original version of our code Janka & 
Hillcbrandt (1989b) performed their analyses for a 
20 M Q progenitor from a core-collapse calculation 
by Hillebrandt (1987). Of course, like our present 
study, these were Monte Carlo simulations on a 
fixed background model, not self-consistent simu- 
lations. Taking into account the different micro- 
physics the mean energies are consistent with our 
present work. The mean energies of v e were some- 
what on the low side relative to v e and the v e lumi- 
nosity was overestimated. Both can be understood 
by the fact that the stellar background contained 
an overly large abundance of neutrons, because 
the model resulted from a post-bounce calculation 
which only included electron neutrino transport. 

Suzuki (1990) studied models with initial tem- 
perature and density profiles typical of proto- 
neutron stars at the beginning of the Kelvin- 
Hclmholtz cooling phase about half a second after 
bounce. He used the relatively stiff nuclear equa- 
tion of state developed by Hillebrandt & Wolff 
(1985). In our table we show the results of the 
model C12. From Suzuki (1991) we took the 
model labeled C20 which includes bremsstrah- 
lung. The model C48 from Suzuki (1993) in- 
cludes multiple-scattering suppression of brems- 
strahlung. Suzuki's models are the only ones from 
the previous literature which go beyond the tradi- 
tional microphysics for transport. It is reassur- 
ing that his ratios of mean energies come closest 
to the ones we find. 

Over the past few years, first results from Boltz- 
mann solvers coupled with hydrodynamic simula- 
tions have become available, for example the un- 
published ones that we used as our Accretion- 
Phase Models I and II. For convenience we in- 
clude them in Table 7. Further, we include a 
very recent accretion phase model of the Garch- 
ing group (Buras et al., personal communication) 
that includes the full set of microphysical in- 
put. Finally, we include two simulations similar 
to the Accretion-Phase Model I, one by Mezza- 
cappa et al. (2001) and the other by Licbcndorfer 
et al. (2001). These latter papers show rms ener- 



gies instead of mean energies. Recalling that the 
former tend to be about 45% larger than the lat- 
ter these results are entirely consistent with our 
Accretion-Phase Models. Moreover, the ratios of 
(e)rms tend to exaggerate the spread between the 
flavor-dependent mean energies because of differ- 
ent amounts of spectral pinching, i.e. different ef- 
fective degeneracy parameters. To illustrate this 
point we take the first two rows from Table 6 as 
an example. The ratio of mean energies for Fermi- 
Dirac spectra with temperatures T\ — 5.2 and 
T 2 = 3.5 and degeneracy parameters rji — 1.1 and 
i]2 = 3.4 is (ei)/(e 2 ) = 1-19, whereas the ratio of 
rms energies equals 1.30. 

To summarize, the frequently assumed exact 
equipartition of the emitted energy among all fla- 
vors appears only in some simulations of the Liv- 
ermore group. We note that the flavor-dependent 
luminosities tend to be quite sensitive to the de- 
tailed atmospheric structure and chemical compo- 
sition. On the other hand, the often-assumed ex- 
treme hierarchy of mean energies was only found 
in the early simulations of Brucnn (1987) and of 
Myra & Burrows (1990), possibly a consequence 
of the neutron-star equation of state used in these 
calculations. 

If we ignore results which appear to be "out- 
liers" , the picture emerging from Table 7 is quite 
consistent with our own findings. For the lumi- 
nosities, typically L„ c L Pa and a factor of 2-3 
between this and L Vfi in either direction, depend- 
ing on the evolutionary phase. For the mean 
energies we read typical ratios in the range of 
(e„ e ) : (e Pe ) : (e„„) = 0.8-0.9 : 1 : 1.1-1.3. The 
more recent simulations involving a Boltzmann 
solvers show a consistent behavior and will in fu- 
ture provide reliable information about neutrino 
fluxes and spectra. 

5. DISCUSSION AND SUMMARY 

We have studied the formation of neutrino spec- 
tra and fluxes in a SN core. Using a Monte 
Carlo code for neutrino transport, we varied the 
microscopic input physics as well as the under- 
lying static proto-neutron star atmosphere. We 
used two background models from self-consistent 
hydrodynamic simulations, and several power-law 
models with varying power-law indices for the den- 
sity and temperature and different values for the 
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electron fraction Y e , taken to be constant. 

The transport opacity is dominated by 
neutral-current scattering on nucleons. In ad- 
dition, there are number-changing processes (nu- 
cleon bremsstrahlung, leptonic pair annihilation) 
and energy-changing processes (nucleon recoil, 
z^e ± scattering). The spectra and fluxes are 
roughly accounted for if one includes one signif- 
icant channel of pair production and one for en- 
ergy exchange in addition to u^N scattering. For 
example, the traditional set of microphysics (iso- 
energetic i/ M N scattering, c + c~ annihilation, and 
v IJ ,e ± scattering) yields comparable spectra and 
fluxes to a calculation where pairs are produced by 
nucleon bremsstrahlung and energy is exchanged 
by nucleon recoil. The overall result is quite ro- 
bust against the detailed choice of microphysics. 

However, in state-of-the-art simulations where 
one aims at a precision better than some 10-20% 
for the fluxes and spectral energies, one needs 
to include bremsstrahlung, leptonic pair annihi- 
lation, neutrino-electron scattering, and energy 
transfer in neutrino- nucleon collisions. Interest- 
ingly, the traditional e + c~ annihilation process is 
always much less important than v c D c annihila- 
tion, a point that we previously raised with our 
collaborators (Buras et al. 2002). None of the re- 
actions studied here can be neglected except per- 
haps the traditional e + c~ annihilation process and 
v^Vc and v^v c scattering. 

The existing treatments of the nuclear-physics 
aspects of the NN — > NNi/D bremsstrahlung pro- 
cess are rather schematic. We find, however, that 
the fluxes and spectra do not depend sensi- 
tively on the exact strength of the bremsstrahlung 
rate. Therefore, while a more adequate treatment 
of bremsstrahlung remains desirable, the final re- 
sults are unlikely to be much affected. 

The transport of and is usually treated 
identically. However, weak-magnetism effects ren- 
der the VfjN and z/^N scattering cross sections 
somewhat different (Horowitz 2002), causing a 
small chemical potential to build up. We find 
that the differences between the average energies 
of and are only a few percent and can thus 
be neglected for most purposes. 

Including all processes works in the direction of 
making the fluxes and spectra of more similar to 
those of v e compared to a calculation with the tra- 



ditional set of input physics. During the accretion 
phase the neutron-star atmosphere is relatively ex- 
panded, i.e. the density and temperature gradi- 
ents are relatively shallow. Our investigation sug- 
gests that during this phase (e„ ) is only slightly 
larger than (e Po ), perhaps by a few percent or 10% 
at most. This result agrees with the first hydro- 
dynamic simulation including all of the relevant 
microphysics except v c v c annihilation (Accretion- 
Phase Model II) provided to us by M. Rampp. For 
the luminosities of the different neutrino species 
one finds Lp e ~ L Uo ~ 2L„ M . The smallness of 
L v is not surprising because the effective radiat- 
ing surface is much smaller than for v c . 

During the Kelvin-Hclmholtz cooling phase the 
neutron-star atmosphere will be more compact, 
the density and temperature gradients will be 
steeper. Therefore, the radiating surfaces for all 
species will become more similar. In this situa- 
tion L v may well become larger than L Po . How- 
ever, the relative luminosities depend sensitively 
on the electron concentration. Therefore, without 
a self-consistent hydrostatic late-time model it is 
difficult to claim this luminosity cross-over with 
confidence. 

The ratio of the spectral energies is most sen- 
sitive to the temperature gradient relative to the 
density gradient. In our power-law models we used 
p cx r~ p and T cx r~ q . Varying q/p between 
0.25 and 0.35 we find that (ep e ) : (e„ ) varies be- 
tween 1 : 1.10 and 1 : 1.22. Noting that the up- 
per range for q/p seems unrealistically large we 
conclude that even at late times the spectral dif- 
ferences should be small; 20% sounds like a safe 
upper limit. We are looking forward to this pre- 
diction being checked in a full-scale self-consistent 
neutron-star evolution model with a Boltzmann 
solver. 

The statements in the previous literature fall 
into two classes. One group of workers, using the 
traditional set of microphysics, found spectral dif- 
ferences between v e and on the 25% level, a 
range which largely agrees with our findings in 
view of the different microphysics. Other papers 
claim ratios as large as (e Po ) : (e„ ) = 1 : 1.8 or 
even exceeding 1:2. We have no explanation for 
these latter results. At least within the framework 
of our simple power-law models we do not under- 
stand which parameter could be reasonably ad- 
justed to reach such extreme spectral differences. 
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In a high-statistics neutrino observation of a fu- 
ture galactic SN one may well be able to discover 
signatures for flavor oscillations. However, when 
studying these questions one has to allow for the 
possibility of very small spectral differences, and 
conversely for the possibility of large flux differ- 
ences. This situation is almost orthogonal to what 
often has been assumed in papers studying possi- 
ble oscillation signatures. A realistic assessment of 
the potential of a future galactic SN to disentan- 
gle different neutrino mixing scenarios should al- 
low for the possibility of very small spectral differ- 
ences among the different flavors of anti-neutrinos. 
The spectral differences between v e and v^, T are 
always much larger, but a large SN neutrino (as 
opposed to anti-neutrino) detector does not exist. 

The diffuse neutrino flux from all past SNe in 
the universe is difficult to detect, although Super- 
Kamiokande has recently established an upper 
limit that touches the upper end of theoretical 
predictions (Malek et al. 2002). If our findings 
are correct, neutrino oscillations will not much en- 
hance the high-energy tail of the spectrum and 
thus will not significantly enhance the event rate. 
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A. MONTE CARLO CODE 

Our Monte Carlo code is based on that devel- 
oped by Janka (1987) where a detailed descrip- 
tion of the numerical aspects can be found. The 
code was first applied to calculations of neutrino 
transport in supcrnovae by Janka & Hillcbrandt 
(1989a,b) and Janka (1991). It uses Monte Carlo 
methods to follow the individual destinies of sam- 
ple neutrinos (particle "packages" with suitably 
attributed weights to represent a number of real 
neutrinos) on their way through the star from the 
moment of creation or inflow to their absorption 



or escape through the inner or outer boundaries. 
The considered stellar background is assumed to 
be spherically symmetric and static, and the sam- 
ple neutrinos are characterized by their weight fac- 
tors and by continuous values of energy, radial po- 
sition and direction of motion, represented by the 
cosine of the angle relative to the radial direction. 
The rates of neutrino interactions with particles of 
the stellar medium can be evaluated by taking into 
account Fermion blocking effects according to the 
local phase-space distributions of neutrinos (Janka 
& Hillcbrandt 1989b). 

As background stellar models we use the ones 
described in Sec. 2.2. They arc defined by ra- 
dial profiles of the density p, temperature T, and 
electron fraction Y e , i.e. the number of electrons 
per baryon. The calculations span the range be- 
tween some inner radius i?j n and outer radius i? out . 
These bound the computational domain which is 
divided into 30 equally spaced radial zones. In 
each zone p, T, and Y e are taken to be constant. 
i?; n is chosen at such high density and temper- 
ature that the neutrinos are in LTE in at least 
the first radial zone. R out is placed in a region 
where the neutrinos essentially stream freely. At 
Ri n neutrinos are injected isotropically according 
to LTE. While a small net flux across the inner 
boundary develops, the neutrinos emerging from 
the star are generated almost exclusively within 
our computational domain. If i?i n is chosen so 
deep that the neutrinos are in LTE, the assumed 
boundary condition for the flux will therefore not 
affect the results. 

The stellar medium is assumed to be in ther- 
modynamic equilibrium with nuclei being com- 
pletely disintegrated into free nuclcons. Based on 
p, T, and Y c we calculate all the required ther- 
modynamic quantities, notably the number densi- 
ties, chemical potentials, and temperatures of pro- 
tons, neutrons, electrons, positrons, and the rel- 
evant neutrinos. The chemical potentials for 
and v T are taken to be zero. Next we compute 
the interaction rates in each radial zone for all in- 
cluded processes. In the simulations discussed in 
the present work, fermion phase-space blocking is 
calculated from the neutrino equilibrium distribu- 
tions instead of the computed phase-space distri- 
butions. This simplification saves a lot of CPU 
time because otherwise the rates have to be re- 
evaluated whenever the distribution of neutrinos 
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has changed after a transport time step. The ap- 
proximation is justified because phase-space block- 
ing is most important in regions where neutrinos 
frequently interact and thus are close to LTE. 

At the start of a Monte Carlo run, 800,000 test 
neutrinos are randomly distributed in the model 
according to the local equilibrium distributions. 
Each test neutrino represents a certain number 
of real neutrinos. In this initial setup the num- 
ber of real neutrinos is determined by LTE. Then 
transport is started. The time step is fixed at 
At = 10~ 7 s; recall that the interaction rates do 
not change. At the beginning of each step neu- 
trino creation takes place. The number of test 
particles that can be created is given by the num- 
ber of neutrinos that were lost through the inner 
and outer boundaries plus those absorbed by the 
medium. Based on At, the production rates, and 
the fact that the inner boundary radiates neutri- 
nos, we calculate the number of neutrinos that are 
produced in one time step and distribute them 
among the available test neutrinos by attributing 
suitable weight factors. The sample particles are 
created within the medium or injected at the inner 
boundary in appropriate proportions. 

During a time step the path of each test parti- 
cle through the stellar atmosphere is followed by 
Monte Carlo sampling. With random numbers we 
decide whether it flies freely or interacts. If it in- 
teracts it can scatter or it can be absorbed; in 
this case we turn to the next particle. For scatter- 
ing we determine the new momentum and position 
and continue with the process until the time step 
is used up. Particles leaving through the lower or 
upper boundaries are eliminated from the trans- 
port. 

After a certain number of time steps (typically 
around 15,000) the neutrino distribution reaches 
a stationary state and further changes occur only 
due to statistical fluctuations. At that stage we 
start averaging the output quantities over the next 
500 time steps. 



B. NEUTRINO PROCESSES 

B.l. Neutrino-Nucleon Scattering 

The rates for the i/ M N reactions are calculated 
following Raffelt (2001). For a neutrino with ini- 
tial energy e\ and final energy e 2 , the differential 



cross section is given by 
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with to — e\ — £2, k the modulus of the momentum 
transfer to the medium, and 9 the scattering angle. 
We do not distinguish between protons and neu- 
trons. Since for nonrelativistic nucleons the scat- 
tering cross section is proportional to C v + iC\, 
the vector current (Cy = — \ for neutrons and 
5 — 2 sin 2 0\y for protons) is small compared to the 
axial component, where we use \Ca\ = 1.26/2. Ne- 
glecting the vector part simplifies the calculations 
significantly and certainly has a smaller effect on 
the scattering rates than other uncertainties, for 
example the in-medium value of the coupling con- 
stants themselves. 

In all of our runs without recoil the structure 
function is given by 



Sno-recoilCw, k ) = 2 7r<5(a>) . 



(B2) 



This corresponds to infinitely heavy nucleons and 
represents the traditional approximation in all pre- 
vious simulations. For the more realistic case of 
recoiling nucleons the structure function becomes 
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with u>k = k 2 /2m. 

Multiplying Eq. (Bl) with the density of nucle- 
ons, ignoring phase space blocking of the essen- 
tially nondegenerate nucleons, yields the differen- 
tial rates that can then be integrated for obtain- 
ing the required energy and angular differential 
rates. In the case of recoil the numerical integra- 
tions are rather tricky because Eq. (Bl) is strongly 
forward peaked. In our code we employ the "re- 
jection method" for obtaining the integrated rates 
(Press et al. 1992). 

B.2. Weak Magnetism 

In order to account for weak magnetism effects, 
we multiply a correction factor to our recoil am- 
plitude, motivated by Horowitz (2002). The cor- 
rection factor has the form 



1± 



4C A (C V + F 2 ) k 
C^(3-cos0) 2m 



(B4) 
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with F2 = 1.019 for protons and F2 = —0.963 for 
neutrons. The upper sign is for neutrinos and the 
lower sign for anti-neutrinos. Expression (B4) is 
always positive; to first order in k/m it is given by 



1± 



4C A (C V + F 2 ) k 
C^(3-cos6>) m 



(B5) 



and thus corresponds to the weak magnetism cor- 
rection of Horowitz (2002), except for one addi- 
tional simplification. In order to keep our pre- 
scription for the recoil-amplitude via the structure 
function (B3) our correction factor must not con- 
tain any energy dependence except for the momen- 
tum transfer k. Therefore, we have substituted 
ei cos 6>, taken in the rest frame of the nuclcon, 
by our momentum transfer k. This is correct for 
forward and backward scattering, but only an ap- 
proximation for other angles. 

B.3. Bremsstrahlung 

We also follow Raffelt (2001). The mfp for 
the absorption of a by inverse bremsstrahlung 
NN is given by 
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The over-barred quantities belong to the that 
is absorbed together with the primary v^. The 
occupation numbers are taken to follow an equi- 
librium distribution with zero chemical potential, 
and I Ca \ = 1-26/2 as in the scattering case. 



B.4. Pair Annihilation 

We now turn to e + e~ 



and , 



The matrix elements for both processes are identi- 
cal up to coupling constants while the phase-space 
integrations only differ by the chemical potentials. 

After summing over all spins and neglecting the 
rest masses, the squared matrix element is 
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(Cv + C A ) 2 u 2 
+ {C v -C A ft 2 \ (B7) 



with the Mandelstam variables t = —2ki ■ k 3 and 
u = —2ki ■ k^. The momenta are assigned to the 



particles as indicated in Fig. 13. The weak inter- 
action constants for c+c~ annihilation are 



-i + 2sin 2 e w , C A = -\ 



while for v e v e annihilation they are 
C\ = Cy = — ■ 



(B8) 
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Fig. 13. — Pair annihilation processes producing 
v^v^ pairs. 

For the interaction rates we have to perform 
the phase-space integrations, using blocking fac- 
tors for the final states and occupation numbers 
for initial-state particles (Hannestad & Madsen 
1995, Yuch & Buchler 1976). Three integrations 
remain that can not be carried out analytically. 

Mu- and tau-leptons are almost absent in proto- 
neutron star atmospheres so that the chemical po- 
tentials of the corresponding neutrinos can be set 
to zero. For the e + e~ reactions the local value of 
= — fi c + = [i c can be obtained from p, T, and 
Y c by inverting 



"c-(Mc) - n c +(fj, c 

^baryons 



(B10) 



where n c -(/i c ) and n e +(fi e ) are Fermi integrals. 
For v c and P e the chemical potential is obtained 
by the relation 



(Bll) 



with the chemical potentials /i p and /i n of protons 
and neutrons, respectively. 

For v c D c annihilation we make use of the fact 
that the energy sphere of lies always deeper 
inside the star than the v c and D e spheres (see 
Fig. 3). Thus v c and v c arc in LTE and are part 
of the medium as far as the transport of is con- 
cerned. This approximation breaks down at larger 
radii where this process is unimportant anyway. 
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For our numerical implementation we normally 
use Fermi-Dirac statistics, but in order to re- 
duce computation time one of the remaining three 
phase-space integrations is approximated by the 
analytic expressions given in Takahashi, El Eid, & 
Hillebrandt (1978). This also requires simplifying 
the blocking factors. With /i c = — /i c + > we can 
approximate the positron occupation number by 
a Maxwell-Boltzmann distribution. For /j, c /T > 2 
this holds to very good accuracy. The greatest de- 
viation is at [i c /T — and yields blocking factors 
too low by about 10%. However, e~ and v e are 
always degenerate in the relevant regions. 

B.5. Scattering on Electrons and Electron 
Neutrinos 

The matrix elements for these reactions are just 
the crossed versions of the leptonic pair processes, 

]T |M| 2 = 8 G 2 F [(C v + C A ) 2 s 2 

+ (C v -C A ) 2 u 2 ] (B12) 



spins 



with the same weak interaction coefficients of 
Eqs. (B8) or (B9) for scattering on e~ or on z/ , 
respectively. For s — 2k\ ■ /c 2 and u = —2k\ ■ ki the 
momenta are assigned to the particles according 
to Fig. 14. Crossing the matrix element Eq. (B12) 
again by interchanging u <-» t, we obtain scatter- 
ing on e + or v c . This is also true for scattering of 
on e~ or i/ c ; scattering of on e + or v c brings 
us back to Eq. (B12). For calculating the rates in 
our code we apply the same approximations as in 
the previous section. 
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Fig. 14. — Leptonic scattering processes. 
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